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Abstract 

The  system  of  hierarchy  equations  for  the  reduced  density 
operators  of  quantiim  statistical  mechanics  is  replaced  by  a 
single  functional  differential  equation  for  a  generating  func- 
tional. A  formal  solution  of  the  initial  value  problem  for  the 
latter  equation  is  obtained^,  leading  to  series  expansions  of 
the  reduced  density  operators.   These  expansions  are  used  to 
obtain  an  improved  derivation  of  the  quantum- mechanical  Boltz- 
mann  equation. 


Table  of  Contents 


Page 


1-   Introduction  1 

2.   Basic  notions  of  quantum  statistical  mechanics  2 

5-   The  Wigner  representation  l4 

k.      The  functional  differential  equation  20 
5-   Derivation  of  the  quantum-mechanical  Boltzmann 

equation  27 

Appendix  51 

Footnotes  and  References  57 


1.   Introduction 

The  statistical-mechanical  treatment  of  a  quantum-mechanical  many- 
particle  system  usually  begins  with  the  'density  operator',  D  ,  which 
is  the  solution  of  an  initial  value  problem  for  an  equation  analogous 
to  the  classical  Liouville  equation.   The  solution  of  this  problem  is 
equivalent  to  the  solution  of  the  quantum-mechanical  n-body  problem  and 
is  therefore  not  practical  when  n  is  large. 

In  order  to  circumvent  this  difficulty  one  may  introduce  'reduced 
density  operators',  F  (s  =  1,2, ...)  which,  in  an  appropriate  representa- 
tion, are  defined  as  integrals  of  the  density  operator.   These  operators 
satisfy  an  infinite  system  of  equation  (in  the  limit  n  -*  <»)  called  the 
'hierarchy  equations',  and  this  system  is,  in  turn,  equivalent  to  a 
single  functional  differential  equation  for  a  'generating  functional', 
L,  which  generates  the  operators  F  . 

In  Section  h   of  this  paper  we  obtain  a  formal  solution  of  the  func- 
tional differential  equation,  and  from  that  solution  we  derive  expansions 
for  the  operators  F  as  power  series  in  the  number  density.   These  expan- 
sions have  the  advantage  that  for  small  densities,  F  may  be  approximated 
by  cL  few  terms  of  the  series.   Then  to  obtain  an  explicit  expression  for 
F  ,  only  certain  k-body  problems,  where  k  is  small,  need  to  be  solved. 
However  the  expansions  probably  suffer  from  the  defect  that  the  remainder 
terms  grow  rapidly  with  time. 

The  work  in  Section  k-   parallels  our  earlier  work  on  the  corresponding 
problem  in  classical  mechanics.   In  fact,  by  introducing  the  appropriate 
representation  of  the  operators,  the  two  treatments  become  formally  almost 


identical.  We  call  this  representation  the  'Wigner  representation' 
because  in  it  the  density  operator  is  represented  by  the  well-known 
'Wigner  function'.   The  Wigner  representation  is  discussed  in  Section  3. 

In  Section  2  we  review  the  basic  notions  of  quantum  statistical 
mechanics  which  we  shall  require.   This  is  necessary  because  these  notions 
are  presented  from  a  point  of  view  which  is  slightly  different  from  cus- 
tomary formulations.   The  unifying  idea  in  this  point  of  view  is  the 
representation  of  Hilbert  space  in  terms  of  a  complete  set  of  commuting 
observables. 

In  Sect: 

Section  k-,    as  the  starting  point  for  a  derivation  of  the  quantum-mechani- 
cal Boltzmann  equation.   By  careful  use  of  operator  representations  we 
are  able  to  avoid  most  of  the  mathematical  difficulties  which  appear  in 
earlier  derivations.   However,  as  we  point  out  in  Section  5,  one  basic 
difficulty  remains  in  that,  at  one  point,  we  are  forced  to  argue  by  ana- 
logy with  the  case  of  classical  mechanics.   The  significance  of  this  dif- 
ficulty and  the  validity  of  the  equation  so  derived  is  discussed  in  the 
concluding  remarks . 

2.    Basic  notions  of  quantum  statistical  mechanics 

Let  V  denote  a  Hilbert  space  of  elements  ^  with  inner  product 
('^,'^  )  and  norm  ||  ^  ||  =  (J,^)^.   The  states  of  a  quantum-mechanical  system 

correspond  to  elements  ^  of  ^  with  |!  ^  ||  =1.   For  many  purposes  it  is 

X  /  N     ^  (1)       (.^) - 

convenient  to  represent  the  elements  £  by  functions  iif(qj  =  i^q   ,...,q 


y'k(q)|%=/k(q)|V^^..dq^^)  <  «..  (2.1) 

This  is  the  well-knovn  Schroedinger^  or  coordinate  representation.  For 
o\ir  purposes  it  is  convenient  to  indicate  the  correspondence  between  ele- 
ments ^  and  functions  \lf(q)  by  writing 

R^d)  =  ♦(q)-  (2.2) 


{q*^',...,,"^') 


If  q  =  <(  q^^',  .  .  .  ;,q^"'  >  denotes  the  complete  set  of  co-ordinate  opera- 


tors then 


R^(q^J^J)  =  q^J^(q);  j  =  l,...,k;  (2.3) 


11111^=  jT  k(q)pciq.  (2.4) 

More  generally,  if  Y  =  <^  Y^  ,  .  .  .  ,Y^        >    is  any  complete  set  of 
commuting  observables  we  shall  denote  the  representation  of  the  Hilbert 
space  in  terms  of  these  operators  by 

V^)  =  c(y)  =  c(y^^\...,j^^^)  (2.5) 

Ry(Y^"^^$)  =  y^'^^c(y);  j  =  l,...,k.  (2.6) 

Here  c(y)  is  a  complex-valued  function.   In  general  ||  ^  ||   will  not  be 
given  by  a  simple  integral  of  the  type  (2.k)   but  may  be  given  by  a  sum 
(discrete  spectrum),  a  combination  of  an  integral  and  a  sum,  etc.   In 
order  to  admit  all  possibilities  we  shall  write 


/|c(y)| 


^din(y).  (2.7) 


With  the  aid  of  the  representation  R^^  functions  of  the  operators 
Y   ,...,Y    are  conveniently  defined  by 

RY[f(y)Jl  =  f(y)c(y)  =  f(y)RYg].  (2.8) 

In  this  paper  we  shall  be  concerned  with  a  system  consisting  of  n 
identical  particles.   Since  an   individual  particle  can  also  be  consider- 
ed as  a  q-uantum  mechanical  system  it  is  essential  to  examine  the  relation- 
ships between  the  Hilbert  spaces  of  the  component  systems  and  the  Hilbert 
space  of  the  composite  system. 

Given  n  identical  quantum-mechanical  systems  with  Hilbert  spaces 
V  ;  V  =  1, . . . ,n:  let  Y  =  Jy^   ,...,Y^  W  be  a  complete  set  of  commu- 
ting  observables  in  X  ,   with  representation 

RY(i)  =  %(y,)  =  c^(y[^^...,y^^));  (2.9) 

V 


Although  the  subscript  v  will  be  used  to  distinguish  the  systems^  we 
choose  the  operators,  and  hence  the  representations,  to  be  the  same  in 
the  n  systems.   Let  X   be  the  space  of  complex-valued  functions 

c(y(^))  =  c(y,,...,yJ  =  c(yW,...,yf))  (2.11) 


for  which 

J     \''(y(n)^\^'^^iyi^-'-M¥j   <".  (2.12) 

Then  X   may  be  viewed  as  representing  a  Hilbert  space  ^   whose  elements 
^  correspond  to  the  functions  c(y/  .) ,    and  for  which 

II  $11^=/  |c(y^^^)!^.bn(y^)...din(y^).  (2.15) 

We  denote  this  correspondence  by 


^      (I)  =  c(y(^^)).  (2.14) 

(n)         ^ 
The  Hilbert  space  M   is  then  said  to  be  the  product   of  the  spaces 


Each  observable,  A  acting  on  ?/      has  a  representation  as  a  func- 
tional operator  (e.g.  a  differential  or  integral  operator)  acting  on 
the  functions  c  (y  ) .   We  now  define  the  operator  A  acting  on  M   by 
requiring  that  the  corresponding  fimctional  operator  acting  on  the 
ftmctions  c(y/  .)   be  the  same  functional  operator,  the  variables  y. 
(j  ^  v)  being  held  fixed. 

The  composite  of  the  n  identical  systems  introduced  above  is  the 
quantum  mechanical  system  whose  Hilbert  space  is  the  product  space  }/ 
and  whose  observables  include  the  operators  A  just  defined,  as  well 
as  all  functions  of  these.   Thus,  in  constructing  the  composite  system 
we  have  identified  every  physical  observable  of  the  component  systems 
with  one  of  the  composite  system.   It  is  easily  seen  that  the  corre- 


spondence  (2.l4)  is  the  representation  of  ^   in  terras  of  the  complete 
set  of  commuting  observables 

Y,  ,  =  /y.....,Y  '\Ax^^\...,1.^''A     .  (2.15) 


(n)  -  ^^l---^n;-|^l 

The  composite  system  defined  ahove  is  called  a  Maxwell-Boltzmann  system 
to  distinguish  it  from  Fermi- Pi rac  and  Bose-Einstein  systems  to  be  dis- 
cussed shortly. 

Let  P  be  any  permutation  of  the  integers  i  1, ...,n  \    .   Then  P 
assigns  to  each  integer  j  in  this  set  another  integer  Pj  in  the  set.   Let 
Py,  s  be  defined  by 

The  permutation  operator  P  on  V  may  now  be  defined  by  means  of  the 
representation  R     : 

R    (P$)  =  c(Py,  J.  (2.17) 

^(n)  ^"^ 

By  means  of  these  permutation  operators^  the  symmetry  operator  J^   and 

the  anti-symmetry  operator  d     are  defined  by 

/„  ^  ^  i:  P   :         «„  =  ^  2  (±)P-         (  =  -^8) 

Here  the  summation  is  over  the  nl  permutations  P,  and  the  positive  or 
negative  sign  is  taken  for  even  or  odd  permutations  respectively.   It 
is  easy  to  show  that  the  permutation  operators  are  unitary  and  that  the 
operators  J     and  (2     are  Hermitian^ 


J   =  J        ,        a   =  1  ,  (2.19) 

and  satisfy  the  condition 

^!=X     ,     <-a^-  (2.20) 


The  last  two  equations  imply  that  J     and  d     are  orthogonal  projectors. 

Let  B  =  B(A^"'"'^)  be  a  fimction  of  all  the  operators  A^"'"'^,  v   =  1,  ...,n; 
i  =  1,  ...,j  :  where  A^"^^  is  an  observable  in  'X'    •   Then  B  is  an  observable 
in  U ,    and  is  said  to  be  a  symmetric  observable  if  B(a)_  ')=   b(A^   )  for 
all  permutations  P  of  Jl,  ...,n|>  .   It  can  easily  be  shown  that  an  opera- 
tor B  is  symmetric  if  and  only  if  it  commutes  with  all  the  permutation 
operators.   From  this  it  follows  that  J     and  d      are  symmetric  and  that 
any  symmetric  operator  B  commutes  with  J     and  d    : 

[^B,JJ  =  0   ;    |b,^J  ^  0.  (2.21) 

Since  J^  (or  ^  )  is  an  orthogonal  projector^  we  may  introduce  the 
subspace  ^z  (or  X,  )  of  )/  onto  which  J  {or  A'    )  projects.   This  is  the 
space  of  all  elements  ^  of  /V  for  which  J  \=^     (or  ^J=  J),  i.e.,  the 
eigen-space  of  J  {ox   IL    )  corresponding  to  the  eigenvalue  1.   The  elements 
^  of  yfr   (or  y^y,    )  are  called  symmetric  vectors  (or  anti -symmetric  vectors) 
and  have  the  property  that  for  every  permutation  operator,  P, 

P^=J     (orP5=+  \).  (2.22) 

Since  every  symmetric  observable,  B,  commutes  with  ^f  ,  it  follows 
that  if  5  is  an  element  of  X'  then  ^\=^J^\=    J  ^\     is  alsc 


30  an 


element  of 


V/  .   Hence  B  is  an  observable  on  }/j    .   The  same  assertion 


is  true  for  -^  .   Thus  we  see  that  there  is  a  natural  correspondence 

between  observables  in  the  component  systems  and  symmetric  f\mctions  of 

these  observables  in  }rj     and  />^  .  A  quantum  mechanical  system  whose 

Hilbert  space  is  f^j    (or  -^  )  and  whose  observables  correspond  in  the 

above  way  to  the  observables  of  the  component  systems  is  called  a  Bose- 

Einstein  system  (or  Fermi -Dirac  system) . 

In  the  remainder  of  this  paper  we  shall  treat  simultaneously  the 

three  types  of  systems  mentioned.   This  is  very  conveniently  doen  by 

introducing  the  operator  Q  ,  where  Q   =  ^  for  Bose-Einstein  systems, 

Q      =    a      for  Fermi-Dirac  systems,  and  Q   =1  for  Maxwell-Boltzmann 
n     n  -^     '       n 

systems.   Thus,  e.g.,   Vq  =  yV  for  the  case  Q   =  1. 

Let  A  =  A(1, ...,j)  be  an  operator  whose  functional  representer  in 
the  representation  R     is  an  integral  operator  acting  on  the  variables 
(l,...,j),  i.e. 

R^    [A(l,...,j)l]  =  /a(y^,...,y.,y',...,y>)c(y-,...,y',y._^^,...,yJ 


(n) 


X  dm(yp...dm(y').  (2.23) 


In  this  case,  we  say  that  'the  operator  A  is  represented  by  the  kernel 

o 
a',  and  indicate  the  correspondence  by  writing 

R,^^^[A(l,...,j)]=a(y^,...,y.,yl,...,yl)=a(y^.),y;j)).  (2.21.) 

In  terms  of  the  notation  introduced  in  (2.23),  A(2, ...,j+l)  is  the 
operator  such  that 

RY^^^[A(2,...,j+l)$]=y^a(y2,...,y.^^,y',...,yl^^)c(y^,y',...,y'_^^,y.^2---^n) 

X  dm(yp...din(y'_^^),  (2.25) 


etc.      It  is  easy  to  see   that  if  A(1, ...,j)    and  B( j+1, . . . ^ j+m)    operate  on 
different   sets   of  variables   then    (a,b]    =  AB  -  BA  =   0,    and 

R  [ab]    =  R  [a]r  [b].  (2.26) 

Let  <  J  y  be  a  complete  orthonormal  set  of  elements  of  V  and  let 
B  be  a  linear  operator  on  V .   The  trace  of  the  operator  B  is  defined  by 

Tr(B)  =  ^  (^^,B  f^),  (2.27) 

V 

provided  the  series  converges  absolutely.   It  can  be  shown  that  the  value 
of  Tr(B)  is  independent  of  the  choice  of  the  set  "S  $  ?  and  if 

%„)M=''(^{n)'='(n)'  f'-^^S' 

then 

Tr(B)  =  J    b(y^^^,y^)din(y^^^)...dm(y^'')).  (2.29) 


It  is  convenient  also  to  define  the  partial  trace  T  ,      .   of  an 

s+1, . . . , J 

operator  with  respect  to  some  of  the  variables.   Thus^  e.g.,  for  s  <  j, 

T  ,      .A(1, ...,j)  is  the  operator  such  that 
s+±, . . . , J 

^^JVl,...,/^l'---^))=  /-(^l---ys'Vl---^J'yi---^^Vl---^j) 

Xdjn(z    )...dm(z_).  (2.30) 


Of  course 


T       A(l, ...,n)  =  TrA(l,  ...,n).  (2.51) 

±,  .  .  .  ,n 


10 


In  quantiim  statistical  mechanics,    the  mean  value   of  an  arbitrary 
observahle,    B,    is   given  by 

<B>    =  Tr(D  B)  (2.52) 

where  D  is  a  non-negative,  Hermitiaxi,  bounded,  completely  continuous 
linear  operator,  called  the  density  operator.   It  is  noiroalized  by  the 
condition 

Tr(D^)  =  1.  (2.33) 

D  varies  with  time  according  to  the  equation 
-itH         itH 


D  (t)  =  e  *    D  (0)  e  ^  (2.5^) 


where  H  is  the  Hamiltonian  operator  of  the  system.   It  follows  from 
(2.5^)  that  D   satisfies  the  differential  equation 

If,  now,  we  have  a  quantum  mechanical  system  with  Hilbert  space 

M     ,    we  must  interpret  the  trace  in  (2.32)  and  (2.53)  as  a  trace, 

Tr    ,  over  X^    ,    i.e.,  the  orthonormal  set  <  ^  (     in  (2.27)  must  con- 

(Q)        ^  ^  ^^ 

sist  of  elements  of  a/^-   It  is  necessary  to  relate  the  trace  Tr,  , 

to  the  trace  Tr  over  M ,    and  this  can  be  done  as  follows: 

Let  <  $  r  be  a  complete  orthonormal  set  in  ^^^^  and  <  If  ?■  a 

complete  orthonormal  set  in  the  ortho- complement  M-^   of  /r     .   Then 

H  H 

the  combined  set  is  a  complete  orthonormal  set  in  /V,  and 

Tr(B)  =  X!  (L^B  L)  +  S  (^,Bf  ).  (2.36) 

V  li     ^^    ^ 


11 


We  assume  that  B  is  an  observable  in  7/  .   (Thus  in  the  cases  0     =    / 
and  Q^  =  A^,    B  is  symmetric.)   It  follows  from  (2.21)  that 

[b,qJ  =  0.  (2.57) 

Now  Q^  $^  =  1^  .  Q^'J'  =  0,  and  Q^  =  Q^.   Thus  from  (2.2?)  and  (2.19) 

Tr(Q^B)  =  Tr(Q^B)  =  Tr(Q^BQ^)  =  ^  t^' V^nt^ 


2   (Qnt'BQ^,  t)=  2   (i^B$^)  =Tr(Q)(B) 


(2.58) 


Thus  in  all  three  cases, 


<B>  =  Tr^^^(D^B)  =  Tr(Q^D^B)  =  Tr(D^Q^B),  (2.59) 

and 


commute  with  Q   .      Therefore 

-itH  itH 


Q^D„(t)    =   e  Vn^°^    *"    "  (2.i^l) 

and 

i^lt   ^Vn^   =    [Wn]-  (2.1.2) 

The  object  of  quantum  statistical  mechanics  is  to  calculate  the  mean 
values  <  B  > .   To  this  end,  one  must  study  the  solutions  of  the  equation 


12 


(2.42).   Although  this  can  be  done  in  abstract  operator  notation,  or  in 
terms  of  any  convenient  representation,  we  shall  work  instead  with  a 
'Wigner  representation'  which  we  shall  introduce  in  the  next  section. 
This  representation  has  the  advantage  of  producing  formulas  which  re- 
semble closely  the  corresponding  formulas  of  classical  statistical  mech- 
anics.. 

The  quantum-mechanical  system  we  shall  consider  consists  of  n  iden- 
tical monatomic  spinless  particles  contained  in  a  finite  volume,  V.   The 
particles  obey  Bose-Einstein,  Fermi-Dirac,  or  Maxwell- Boltzmann  's.tatis- 
tics',  according  as  0  =  J  ,    Ql     =    (2   ,    or  Q     =   1.      The  Hamiltonian  of 
the  system  is  given  by 

n 


"    P"     ^  ^  (P?^' 


T  = 


s  ;i^ss^. 


.  -,    2m     ^  ^   2m 
x=l  1=1  o^l 


U  =       2      <t>(|li-qJ)-  (2-^5) 

1  S  i  <  j  S  n     ^   ^ 

Here  m  denotes  the  mass  of  a  particle,  1^   =   <,   <1^}1^>(1^  J   aj^d 

p.  =  J  p.,p.,p?  f   denote  the  cartesian  coordinates  and  the  conjugate 

moments  of  the  i   particle,  (j)  is  the  inter-particle  potential,  and 

u^  (q.  )  is  the  potential  due  to  the  containing  boundary.   Thus  Uy(q.)  is 

constant  inside  V  and  rapidly  approaches  infinity  at  the  boundary.   At 

an  appropriate  point  we  shall  let  n  and  V  tend  to  infinity  in  such  a  way 


15 


that  V  =  v/n  approaches  a  finite  limit  (the  specific  volume).   In  this 
process,  the  effect  of  the  terms  u^,(q.  )  vanishes   ;  hence  we  shall  omit 
these  terms  from  further  equations.   The  resulting  Harailtonian  is 

(2.46) 

In  our  analysis  we  shall  make  use  of  Hamiltonians 

T^;   s  =  1,2,...  (2.J+7) 

corresponding  to  systems  of  s  interacting  and  non-interacting  particles. 

In  what  follows  we  shall  make  frequent  use  of  a  representation  in 
terms  of  the  cojnplete  set  of  moment\im  operators 

Since  this  representation  will  be  used  so  frequently  we  rfill  denote  it 

simply  by  R,  rather  thsui  R    .   The  momentum  representation  is  charac- 

^(n) 
terized  by 

R(l)  =  <J'(P(n))'  (2.49) 

R(P^$)  =  v^HvJ,  (2.50) 

\\lf=   /l4)(P(,))  !%(„)•  (2-51) 

Here  dp,  ^  =  dp  ...dp  =  dp  ...dp  ,  and  the  integration  is  over  all  values 
°^P(n)- 


Ih 


3.    The  Wigner  representation 

Let  A(1^...,s)  be  an  operator  such  that 

r[a(1,...,s)]  =  a(p^^),pj^j).  (3.1) 

Here,  as  in  Section  2,    R  denotes  the  momentum-representation.   Let 


^^^(s)^p(s)^  =  ^^"^^'^'Z  "^^r-fSs)-  "(s) 


X  ^^P(s)-  "(s)'P(s)-'  '(s)^^^(s)-  ^5.2) 


It  follows,  by  inversion  of  the  Fourier  transform,  that 


-(P(s)^P(s))  =  j   "^   I^(s)-^P(s)-P(s)) 


■[^(s)^  i^p(s)^p;s))]^^(s)-  ^5.5) 


These  equations  define  a  correspondence  between  functions  a(q,  wP/  \)  and 
operators  A(1, ...,s).   It  is  not  apparent  that  a  is  the  representer  of  the 
operator  A(1, ...,s)  in  terms  of  any  complete  set  of  commuting  observables. 
Nevertheless  it  is  convenient  to  call  a  'the  representer  of  A  in  a  Wigner 
representation '  and  to  indicate  the  correspondence  by  writing 

w[a(1,...,s)]  =  a(q^^^,p^^^).  (3.^) 

With  the  aid  of  (2.26)  it  is  easy  to  see  that 

w[a(1,  .  . .,  j)B(  j+1,  .  .  .,  j+mQ  =  w[a(1,  .  .  .,  j)]w[B(  j+l,  .  .  .,  j+m)]  . 

(3.5) 
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For  convenience,  we  shall  sometimes  denote  the  pair  q,p  by  x.   Thus,  e.g., 

^(s)'P(s)  =  ^(s)>   V..-,q3,Pi,.-.,P3  =  x^,...,x^;      (5-6) 

etc. 

In  terms  of  the  above  definitions,  the  well-known  Wigner  function, 
w  (t,X/  v)  is  simply  the  Wigner- representer  of  the  operator  Q  D  (t): 

V^'-(n))  =  V^'^(n)^P(n)^  =  ^[Vn^^^]'  ^^'^^ 

The  properties  of  the  Wigner  function  are  discussed  extensively  in  [llj  . 
Here  we  shall  outline  briefly  those  which  we  shall  require. 

Since  Q  and  D  are  Hermitian  and  commute,  it  follows  that  their  pro- 
duct is  Hermitian  and  hence  that  w  is  real,  although  not  necessarily  non- 
negative.   Given  an  arbitrary  function  g(q/  -.,j>,    -.) ,    it  is  well  known 
that  because  of  ambiguities  arising  from  the  ordering  of  non-commuting 
operators  there  may  exist  more  than  one  operator  corresponding  to  g.   How- 
ever, among  these  operators  one  may  specify  a  unique  operator,  called  the 
Weyl  assignment.   If  G  is  the  operator  corresponding,  under  this  assign- 
ment, to  g(q/  \,P/  \)  then 
(n)   (n) 

<G>=  T^(Q^D^G)  =y^wJt,q(^^,P(^pg(q^^^,p^^^)dq^^^dp^^^. 

(5.8) 

The  Weyl  assignment  and  the  proof  of  (5-8)  are  given  in  [llj  • 

In  classical  statistical  mechanics,  the  mean  value  of  a  phase 
function  g(q/  x,^,    -, )  is  given  by  a  formula  identical  to  (3'8)  where 
w  is  the  phase  space  probability  density  function,  and  satisfies  the 
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Liouvllle  equation  (see  [l] ) .   The  great  virtue  of  the  Wigner  function 
formalism  Is  this  similarity  to  the  classical  case.   Indeed  the  Wigner 
function  satisfies  a  'quantum-mechanical  Liouvllle  equation',  similar 
to  the  classical  equation  (but  differing  in  the  crucial  fact  that  it  is 
not  of  first  order) . 

This  equation  is  simply  the  Wigner  representation  of  (2^2)  and  is 
given  by 


^    I  "  "/     1  s  t^k  g  n 

{w  ;T  1 
\    n'    nj 


.,  'w  =  0. 
jk  n 


(5.9) 


Here  ■^  w  ;T_  >  denotes  the  Polsson  bracket,  and  the  operator  ©  .,  •   is 
given  by 


9  .  -w  =  -  I  sin  I  V 
Jk  n     ^         \2  P 


(„,  \„,j«IV%i'"n(*''-(n)'%)> 


{3-10) 


We  note  that 


-V    (t)(k.-q,  |)-V    w  (t.q,  wP/  .)    +   0(^  ) 
^(n)    -^  ^    P(n)    '^'^^'^(^^' 


=  -  {(Kl^j-lkD'^n)  ■^°^^^)' 


(3.11) 


The  classical  Liouvllle  equation  is  given  by  (3-9)  with  ©  •j^-'^n  ^'^placed 
l^y  -  S 't'C  K-'^lv  I  ^'^  I  '    ^^^^  ■'■*  ^^    clearly  the  limiting  form  of  the 
quantum-mechanical  equation  for -h  -*  0.  Various  forms  of  the  operator 

6   .   are  given  in  llll . 
jk       °        ^    -' 


17 


With  the  aid  of  the  Wigner  representation,  ve  can  introduce  a  formal 

solution  of  (3'9)  derived  from  (2.^1).   This  solution  is  given  in  terms 

(k) 
of  a  solution  operator  S(l, ...,k)  defined  as  follows:   If 

t 
a(q/  wP/  \ )  =  w[a(1,  .  .  .,s)J,  and  k  S  s  then 


s('')(l,...,k)a(x(^,)  =  s('^'(l,...,k)a(l(^,,p,^)) 


e  A(l, ...,s)e 


(5.12) 

Here  H  is  defined  by  (2.k'~[) .      The  solution  operator  s|  ^°'^(l,  .  .  .  ,k)  is  de- 
fined by  (3-12)  with  R     replaced  by  H    .   In  an  analogous  way  we  can  define 

(k) 
the  operator  S^   (2, ...,k+l),  etc.  Where  no  confusion  can  arise,  we  omit 

the  explicit  indication  of  the  variables.   Thus,  e.g.. 


4'^^('l(s)^P(s))=  s['^(l-'-^)-('l(s)'P(s)^ 


(3.15) 


It  follows  immediately  from  the  definitions  that  S 


(k). 


(k)   (k)     (k) 


^t   ^t 
\        ^2 


(5.1M 


From  (2.4l)  it  is  clear  that  the  formal  solution  to  the  initial  value 
problem  for  the  quantum-mechanical  Liousvile  equation  is  given  by 


f   {t,qr    N,p^  x)  =  S*^^^  w  (0,q/  .,p,    J. 
n^  '^{n)^[n}'  -t   n^  -"^(n)  -^(n)' 


(5.15) 
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For  use  in  Section  5^  it  will  be  convenient  to  define  an  operator 
2^  which,  in  the  Wigner  representation,  corresponds  to  multiplication 
of  an  operator  by  the  operator  Ql.      Thus  if  a(q,  s,p,  .)  =  w[a(1,  . .  . ,  s)^ , 
let 

^k"(^(s)^P(s))="[V(l---^^]-  (5.16) 

In  the  case  of  classical  statistical  mechanics  solution  operators 
may  also  be  defined.   One  first  introduces  the  solution  operator  for 
the  system  of  Hamilton's  equations  corresponding  to  the  Hamiltonian  K 
(or  R   ).   If  we  denote  the  solution  of  these  equations  by 
X,  ,(t)  =  <  q/,  N(t),p/  s (t)  y    ,    then  this  solution  operator  is  defined  by 

4'*  -(.)(«)  =  ''«(')■  (5-^^) 

This  operator  maps  each  point  in  6k-dimensional  phase  space  into  its 
image  after  time  t  along  the  classical  trajectory.   The  classical  solu- 
tion operator  for  functions  is  then  defined  by 

sf)(l,...A)a(x(^^)  =  a[s(^)x^^),x^^^,...,xj,  (5-18) 

sf)(2,...,k+l)a(x(^))  =  a[x^,sf){x2,...,x^^^),x^^2----s]'   (5 
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etc.   If,  now,  we  interpret  S^^^    as  the  classical  solution  operator,  the 
solution  of  the  initial  value  problem  for  the  classical  Liouville  equation 
is  given  by  (3-15) • 
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It  is  important  to  observe  that  in  the  case  of  free  particles  (())  =  0; 
H  =  H   )  the  classical  and  quantum-mechajiical  solution  operators  coincide. 
This  is  seen  most  easily  from  the  fact  that  in  this  case  the  two  Liouville 
equations  are  identical. 

In  both  the  classical  and  the  quantum-mechanical  cases,  the  determina- 
tion of  the  operator  s|   is  eqiiivalent  to  the  solution  of  the  corresponding 
mechanical  problem  and  hence  is,  in  general,  not  practical  for  n  >  2.   In  an 
attempt  to  circumvent  this  difficulty,  we  introduce  's -particle  Wigner  func- 
tions '  defined  by 


f   (t,x,  J  =  f   (t,x,,...,x  )  =  V^  /  w  (t,x,,...,x  )dx  ,...±jc  ; 
n,s^  '    (s)'    n,s^  '    ±'        '    s'  J        n^  '  1'    '    n'      s+1     n' 

s  =  l,2,...,n.  (5-20) 


By  integrating  equation  (3-9)  with  respect  to  x   ,...,x  ,  one  obtains  the 


af 

n,  s 


(f   ;T  }  +  ^,      e.'f 
{_   n,s'  sj    -,^'TT'^   Jk  n,s 


•f    ,(t,x,,...,x  Jdx  ,j 
s+1  n,s+l^    1'     s+1'   s+1' 


s  =  1,2, ...,n-l.  (5.21) 

Here  T  is  defined  by  (2.4ij-)  and  v  =  v/n.  We  observe  that  the  left  side  of 
(5.21)  is  identical  to  (5-9)  for  s  =  n.   In  the  derivation  of  (5.21)  one 


makes  use  of  the  symmetry  of  w  ,   which  is  a  consequence  of  the  symmetry 
of  the  operator  Q  D  .   The  details  of  the  derivation  are  given  in  \ll\ . 

In  the  next  section  we  shall  obtain  a  formal  solution  of  the  initial 
value  problem  for  the  system  of  equations  (5. 21)  in  the  limit  V  ^  oo^ 
n  -♦  00^  v/n  =  V.   In  the  following  section  we  shall  use  that  solution  to 
derive  a  quantum-mechanical  analogue  of  the  Boltzmann  equation. 

k.        The  functional  differential  equation 

The  work  in  this  section  parallels  the  development  of  the  classical 
case  given  in  Sections  2  and  3  of  [l] .   In  fact,  if  we  simply  replace  the 
classical  density  function  D  of  [l]  by  the  Wigner  function  w^,  and  the 
operation  j  ({)(  |  q  .-q,  | )  J  \      hy  -  G   •  ,    the  analysis  becomes  identical. 
Hence  we  shall  only  outline  the  development  here,  and  refer  to  [l]  for  the 
details . 

We  begin  by  introducing  the  functional 

n 
Ljt,u]  =jr,jt,x_,...,xj  n_^   [l  +  w(x.)]dx^...dx^,         (^.1) 

which  is  defined  on  the  domain  of  functions  u(x)  for  which  the  integral 
converges .   By  functional  differentiation  we  obtain 


6u(x^)...Su(Xg) 


U;=0 


n  (n-s) 


E"   (t,x,,...,x  );   s  =  0,1.2,. ..,n.     {h.2) 
n,  s^  ^  1^    '    s" 


With  the  aid  of  (^^-.2),  L  may  now  be  expressed  as  a  (finite)  series  expan- 
sion around  u  =  0. 


[t,u]  =1+    S    Tr(l-^)---(l  -^) 
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Ix  /.  S-1, 


/  \,it,^^>.-.,^^M.^)..M^^H^^...^^.  (^.5) 


A  differential  equation  for  L     may  be  obtained  by  multiplying   (3-9) 
n 
by     IT    [l  +  vu(x.)J    and  integrating  with  respect   to  x   , ...^x   .      We  then 

i=l  ^  in 

let  n  -►  00  ^   V  ^  00     in  the  resulting  equation.      In  this  process  we  always 

assume  that  v/n  =  v  approaches   a  finite   limit.      Thus    we   obtain 


/[u(x,).^][T(P,);^]dx^ 


Here 

2 


T(p)  =  ^  and  LQt,u]    =      lim     L^[t,u]  .  (1^.5) 

n  ->  00 
V  ^   oo 

lim 

If  we   set  f  (t.x,  ,  . .  .,x   )   =  n^oo   f        (t.x^ ,  . .  .,x   )   and  pass   to  the  limit 
s^    '    1-^         -^    s'  _^     n,  s^    '    1'         '    s'  ^ 

in   (5-21)    and   (4. J),   we   obtain 


-      ^         /    ©•  n-  f      -,("t;X,  ,  .  ..,X      Jc 

V      -^^        /       1,S+1      s+1^    '    1'  '    s+1' 

1=1     ^        ' 


s  =  1,2,...  (k.6) 
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L[t,u]  =  1  +  S  jT   f    f3(t,x^,...,x^)u(x^)...u(xjdx^...d^^. 


(^•7) 


It  follows  from  (4.7)  that 


f  (t,x^,...,x  )  =  —J A^V-7 7 

s'  '  1^    ^  s^    5u(x-,  ) . .  .5u(x  ) 


(^.8) 


The  equations  of  the  infinite  system  {h .6)  are  called  the  quantum- 
mechanical  hierarchy  equations.  The  functional  L  is  called  the  genera- 
ting functional  because  it  generates  the  functions  f  by  means  of  (4.8). 

In  order  to  solve  the  initial  value  problem  for  (4.4),  we  begin 
by  examining  the  case  of  zero  density,  l/v  =0.  We  shall  use  superscripts 
'o'  to  denote  this  case.   It  is  convenient  now  to  denote  the  independent 
function  by  w  instead  of  u.   Thus  from  -(4.7),  (4.4)  and  (4.8), 


'[t,w]  -   1  +     ^    ~T    J   f°(t,x^,...,x^)w(x^)...w(x^)da^... 


...dx  , 
s 

(^.9) 


/v(x^) 


T(pJ; 


1^'  5w(x^) 


i  /"("i)"^"2^\2-5w(x'MxJ  ^1^2  =  ^^  (^-lO) 


25 


fr(t,x,,...,x J 


!'•••'    s^~      6w(x^)...6v(Xg) 


(^.11) 


and  {^.6)   reduces  to 


s 


■>J  ^  ,,s 


lSj<kSs 


jk     s 


1,2,., 


(^.12) 


Comparing  (^.12)  with  {3 •9)   we  observe  that  f  satisfies  the  quantum- 
mechanical  Liouville  equation  for  s  particles.   It  follows  from  (3. 15) 
that  the  solution  of  (4.12)  is  given  by 


f°(t,x^,...,x^)=s(^)  f°(0,x^,...,x^). 


1.2, 


(^.15) 


where  s\_"     is  the  quantum-mechanical  solution  operator  for  a  system  of 
s  particles.   On  inserting  (4.15)  in  (4.9)  we  obtain  the  solution  of 
(4.10)  subject  to  the  initial  conditions 

L°[o,w]  =  1  +  X!  FT  /   f°(0,x^,  ...,x^)w(x^)...w(x^)dx^...d^g. 
s=l    ■  ^ 

(^.IM 


The   f  (0,x  , ...,x  )  are  the  given  initial  data. 

In  order  to  solve  the  general  equation  (4.4)  we  compare  it  with 
(4.10).   This  comparison  suggests  that  we  try  a  solution  of  the  form 


.[t,uj  =  L°[t,w];   w(x)  =  u(x)  +  -. 


(^.15) 


2k 


Inserting  this  expression  in  {k .k)   we  see  at  once  that  the  equation  is 
satisfied  by  virtue  of  the  fact  that  L  satisfies  (ij-.lO). 

But  (^.'+)  must  be  solved  subject  to  the  initial  conditions 

l[o,u]  =  ^[u]  =  1+  ^   ^  /  f  (0,x  ,  ...,Xg)u(x^)...u(Xg)dx^...dXg. 
s=  1   '    "^ 

(^.16) 

The  functions  f  (0,x  , ...,x  )  are  the  given  initial  data.   In  terms  of 
(4.15)  the  initial  conditions  become 

L°[o,v]  =  ^[u];   u  =  V  -  i.  (4.17) 

Now  let  L°[t,w]  be  the  solution  of  the  initial  value  problem  for 
(4.10)  with  initial  conditions  presented  in  (4.17).   Then  LJjtjuJ,  defined 
by  (4.15),  is  the  solution  of  the  initial  value  problem  for  (4.4)  with 
initial  conditions  presented  in  (4.l6). 

One  may  now  obtain  expansions  of  the  fimctions  f^(t,X-|^,  .  .  .  ,x^)  in 
powers  of  the  density,  l/v,  as  follows:   From  (4.8)  and  (4.15), 


6\[t,u] 
f^(t,x^,...,x^)    =   5^(^^)...5^(^J 


5V[t,w] 
5w(x^)...5w(Xg) 


(^.18) 


If  we  apply  5^/Sw(x  ) . . .5w(x  )  to  (4.9)  and  set  w  =  l/v,  we  obtain 

f  (t,x  ,...,x  )  as  a  power  series  in  l/v  with  coefficients  involving  the 

functions  f°(t,x^, . . .,x^)  =  S^^^  f°(0,x^, . . . ,x^ ) .   Now  from  (4.11)  and 

(^.17) 
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f3(0,x,,...,x^) 


6V[o,v] 
5w(x^)...6w(x^) 


5V[u]  I 

Su(x^)...Su(x^)  ^ 


(^.19) 


If  We  apply  5  /6u(x-j^) .  .  .6u(x^)  to  (^.l6)  and  set  u  =  -l/v,  we  obtain 

(0,x  , ...,x  )  as  a  power  series  in  -l/v  with  coefficients  involving 
the  fimctions  f^(0,x  ,  ...,x  ),  which  are  the  given  initial  data.   If  we 
now  insert  the  series  expression  for  f°(0,x-^ ,  .  .  .  ,x  )  into  the  series  ex- 


pression for  f  (t,x^...,x  )  we  obtain 


f  (t,x,, ...,x  ) 
s^    1    '  s' 


k=o      ^ 


^  j:(k-j):  s-t  ^^' 


j+s) 


xs(^-J'°)(j..s+l,...,lc^s)f^_^^(0,x^,...,x^_^^) 


1,2, 


(4.20) 


The  series  expansion  (^+.20)  for  f  has  the  advantage  that  for  small 
densities  (-  «  1)  the  fimction  f  is  approximated  by  terminating  the 
series  after  a  few  terms.   Now  the  functions  f  of  main  interest  are  those 
for  which  s  is  small,  and  for  these  functions,  the  calculation  of  the  lead- 
ing terms  of  the  expansion  reqtiires  a  knowledge  only  of  solution  operators 

(k) 
S)        where  k  is  small,  and  initial  data  f . (0,x  , . . .,x . )  where  j  is  small. 
"^  J    1      J 

The  expansion  (4.20)  probably  suffers,  however,  from  the  same  defect  as 
the  corresponding  expansion  in  the  classical  case,  namely  that  the  remain- 
der terms  grow  rapidly  with  time. 
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In  order  to  obtain  the  operator  analogue  of  (^.20)  we  must  first 
introduce  's-particle  density  operators' 


^   (t,l,...,s)  =  V^  T  ^,       [QDCt)],  (^.21) 

n, s^  ^  ^  s+l^...,n|_n  n^  'J'  ^  ' 


and 


F  (t,l, ...,s)  =  lira   F   (t,l,...,s) 
n  -♦  00 


(4.22) 


Then  it  is  easy  to  show  that 

,  f^(t,x^....,x^)  =w[Fjt,l,...,s)] 
and  the  analogue  of  (4.20)  is 


(4.23) 


F  (t,l,.-.,s) 


t  (|>^ 


s+1, .... s+k 


V  (-1) 
Zj  TT7k 


k-J 


j;(k-j)' 


xexp^-  ^[h._^Ji,...,j+s)  + 


H^_.^Jj+sH-l,...,k+s)]l 


xF^^JO,l,...,k+s)exp^  f[H._^^(l,...,j-»-s) 


+  \_.^^{j+s+l,...,]^+s)j> 


s  =  1,2,  ...  .       (4.24) 
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In  the  case  of  classical  mechanics,  s-particle  density  functions 
are  defined  as  in  (5-20),  where  w  is  replaced  by  the  classical  density 
fimction.   They  satisfy  a  system  of  hierarchy  equations,  and  the  solu- 
tion of  the  initial  value  problem  for  these  equations  is  given  by  ('+.20) 
where  the  solution  operators  are  those  appropriate  to  classical  mechan- 
ics.  These  matters  are  discussed  in  detail  in  [l] .   In  that  paper  we 
applied  the  leading  terms  of  the  expansion- (^.20 )  for  s  =  1  to  obtain 
a  simple  derivation  of  the  Boltzmann  equation.   In  the  next  section,  we 
shall  attempt  to  derive  the  quantum-mechanical  analogue  of  the  Boltzmann 
equation  in  a  similar  manner. 

5.    Derivation  of  the  quantum-mechanical  Boltzmann  equation 

1^ 
In  1955;  Uehling  and  Uhlenbeck   suggested  that  certain  modifications 

of  the  collision  integral  of  the  classical  Boltzmann  equation  would  be 

appropriate  to  the  case  of  quantum  mechanics.   Since  that  time  various 

,,       10,15,16,17,18,19,  +^      ^   -,   ^  -^    +V,        T^--.       ^     ^V, 

authors  '    ^'  have  attempted  to  verify  the  validity  of  the 

resulting  equation  by  deriving  it  from  the  basic  formulas  of  quaxitum  sta- 
tistical mechanics.   In  our  opinion,  none  of  these  derivations,  including 
the  one  to  be  presented  here,  is  entirely  satisfactory.   Some  of  the  in- 
consistencies of  the  earlier  derivations  can  be  removed  by  proper  use  of 
operator  representations  and  application  of  recent  developments  in  time- 
dependent  scattering  theory.   The  purpose  of  the  present  derivation  is  to 
remove  these  inconsistencies  and  to  isolate  what  we  consider  to  be  the 
real  difficulty. 

At  one  point  in  our  derivation  we,  like  the  earlier  authors,  are 
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forced  to  argue  by  analogy  with  the  classical  case.   It  is  at  this  point 
that  we  feel  the  basic  difficulty  remains.   After  the  derivation  is  com- 
pleted we  shall  add  some  concluding  remarks  about  this  difficulty  and  the 
possibility  of  removing  it. 

The  starting  point  of  our  derivation  is  the  solution  (i|.20)  of  the 

quantum-mechanical  hierarchy  equations.   Following  the  suggestion  of 

20 
Kirkwood   and  others,  we  shall  make  use  of  a  time  average  of  f,  .   The 

form  of  the  'molecular  chaos  assumption'  appropriate  to  each  of  our  three 

cases  (Q  =J,Q=<!Z,Q=1)   may  be  stated  in  terms  of  density 

operators  as  follows: 

FgCtA^S)  =  Q2F^(t,l)F^(t,2).  (5.1) 

The  justification  of  this  assumption  is  discussed  in  [lo] .  We  shall 
assume  that  it  holds  for  all  t.   It  follows  from  (3-5);  (3-l6),  and  (^.23) 
that 

fgCt^x^.Xg)  =  ^^f^{t,x^)f^{t,^^).  (5.2) 

Thus  (5.2)  is  the  Wigner  representation  of  the  molecular  chaos  assumption. 
If  we  set  s  =  1  in  (i)-.20),  we  obtain 

f^(t,x^)  =  s^j)f^(0,x^)  f  i  /[s^^)f2(0,x^,X2)-S^j)(l)S^i\2)f2(0,x^,X2)]dx. 


(5.3) 


Since  initial  values  may  be  specified  at  an  arbitrary  time  (say  at  t  -  2) 
we  obtain J  by  replacing  t  by  t. 
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r^(t  .  1  x^)  =  s(^)f^(t  -  l,x^)  .  i  f  \s[ff^it  -  l,x^,x^) 

-   S^^V)S^^^(2)f2(t   -  |.x^,X2)]dx2  +  of  i^^   .  (^.k) 


In   a  similar  manner  we   obtain 

f^(t   -  l,x^)   =  s[^)f^(t,x^)    +  0(i);      S^^)f^(t+s,x^)   =  f^(t,x^)   +  0(i) . 
2 

(5.5) 

If  we  apply  the  operator  S^  "^(l)  to  (5-^)  the  result  is 
2 


;[^)f^(t  +  i,x^)  -  S^^^f^(t  -  i  x^) 


(5-6) 


i   r  S^^)(l)S^ff2(t  -  l,x^,X2)-S^^^(l)S^^)(2)f2(t-|,x  ,X2)dx2+o[^ 


(5-7) 


We  now  introduce  the  transformation  of  integration  variables, 


x_  =  S^]l   x'  in  (5.7)-   Since  s;"^^x  =  si'^^Cq^p)  =  (q  +  -  p,p)  we  see  that 
dx^  =  dXp.   If  we  then  change  x'  to  x^  in  the  resulting  integral  and  note 


2     2 
that 


2  "   2 

we  obtain 


i  I  S^2,o)3(2)^^(^_x^^^^^^)_3(l)^^)^(l)^2)f2(t-i,x^,X2)  ^^g-^ofi^j 


(5.8) 
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From  the  molecular  chaos  assumption  (5.2),  and  from  (5>5) 

=^2S^l)(l)f^(t,x^)sJ^)(2)f^(t,X2)  +  0(i) 


.?2S^^'°^fl(t,x^)f^(t,X2)  +  0(i), 


(5.9) 


Thus,  since  5p  commutes  with  s)_        and  s|  "^  ' 


Z  = 


,(2,0)3(2)^(2,0)^  ^^^^^^^^)^^^^^^^)_^^^^(^^^_^)^^^^^^^) 


(5-10) 


We  now  introduce   a  time   average,    f,    of  the   one-particle  Wigner  fiinction, 


f(t,x^)    =  i      J      S^^^f^(t+s,x^)ds, 


(5.11) 


and  we  define  the  total  time  derivative  of  a  function  h(t,x  )j 


V^*>^i)  =  h  si^^h(t+z,x^; 


(5.12) 


D^h(t,x^)  =_h(t+z,q^+^p^,p^) 


[Ie-^i-\]^-    (5-^5) 
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From  (5-11),  (5.12)  and  (5.6), 
2 


^r^v^^i-i)-^^i^M-J-i) 


(5.1^) 


and  from  (5-11)  and  (5.5) 

f^(t,x^)  =  f(t,x^)  +  0(i).  (5.15) 

If  we  insert  (5-15)  into  (5. 10),  and  the  latter  in  (5-1^)  we  obtain 


D,f(t,x^)  =  i/ 


s['''°h[fs[^'°>i^f(z,^^)t{t,^^)-!^r(t,y.^)t(t,^^) 


^H;^ 


dx. 


(5.16) 


We  recognize  the  left  side  of  (5 -16)  as  that  of  the  Boltzmann  equa- 
tion.  Our  object  now  will  be  to  transform  the  right  side  into  the  Uehling- 
Uhlenbeck  modification  of  the  collision  integral. 

So  far  the  analysis  in  this  section  is  valid  for  the  case  of  classical 
mechanics  as  well  as  for  the  case  of  quantum  mechanics,  provided  we  inter- 
pret the  solution  operators  according  to  their  classical  definition,  and 
take  <2p  =  1.   Let  us,  for  the  moment,  assiome  that  we  are  dealing  with  the 
classical  case  and  examine  the  solution  operators  appearing  in  (5-l6).   The 
following  discussion  depends  on  the  fact  that  the  classical  solution  opera- 
tor is  defined  first  for  points  in  phase  space  (j-l?)  and  then  for 
functions  (3.I8).  With  the  aid  of  Figure  1  we  examine  the  effect  of  the 
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^ 2^|PZ-P,I 


Fiflure  |:    CONFIGURATION   SPACE   FOR  CLASSICAL   BINARY 
COLLISION. 
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operator  ^\ /'^     ^--i  ^\ /'^       ""^  Points  ix^^Xgi  =  [  <li^Pi^  l2'P2  j  ^°^ 
fixed  values  of  P-,>Pp,  and  q  ,  and  various  choices  of  qp. 

In  Figure  1,   A,B,C  together  form  a  cylindrical  region  with  spheri- 
cal end  surfaces.   The  radius  of  the  cross  section  of  the  cylinder  is  r  , 
the  interaction  distance  of  the  interparticle  potential;  i.e.,  we  assume 
that 

(5. IT) 


,(2,o; 


For  points  q^  in  C,  the  operator  S^ /p   maps  qp  to  q',  the  position  of 


S^  /p  qp  in  a  coordinate  system  in  which  q  is  fixed.   Under  S  /p   ,  ^p 
moves  relative  to  q  ,  a  distance  p-  [pp-P-,  |  in  the  direction  p  -  p  , 

along  a  trajectory  which  is  a  straight  line  because  S  /p   corresponds  to 

(2) 
a  Hamiltonian  Hp   with  no  interaction.   Then  S    maps  q'  to  qp  along  a 

trajectory  which  is  first  straight,  then  curved  in  the  interaction  region 
(the  sphere  of  radius  r  about  q  )  and  then  straight,  with  direction 
-(Pp~  Pt ) •   Finally  S  /p   maps  q^  back  along  tne  straight  trajectory  with 
direction  p  -  p  to  q'"  .   The  momenta  P-|;Pp  are  unaffected  by  the  opera- 
tors S  /p   but  are  transformed  into  P-|j>Pp  hy  S^   .   Now  if  t  is  increased, 
the  points  q'  and  q"  move  outward  along  the  straight  trajectories  but  q'" 
does  not  change  its  position.   Thus  for  points  q  in  C, 


(^.°)s(fs(^'°[x^,x,j=  li„^sf'°)s.(  =  )s(^'°){x,,xj.  (5.18) 


It  is  easy  to  see  that  points  outside  the  region  A,B,C  simply  map 
into  themselves  under  the  operator  S  /^  S  S  /I  .  (Of  course,  the 
length  of  the  cylindrical  region  increases  linearly  with  t.)   It  follows 


5^ 


that  the  integrand  in  (5.I6)  is  zero  for  points  q  outside  the  regie 
A,B,C,  hence  certainly  for  points  in  the  region 

k-p|  >i^  |P2-  pj  +  ^1 


where 


i2   ^1 


^2   ^1 
IPg  -  Pj 


Therefore,    we    can  replace   f(t,x   )f(t,x    )    in    (5.I6)   by  a  fimctic 

g('t,q-[_,q2'Pi^P2^   where 

f  f  for    |q.p|    S  ^    Ip^-P^I    +  r^ 


0      for    |q-p|    >^    IP2-P1I    +   ^1 


(5.19) 


(5.20) 


(5.21) 


f(t,q^,p^)  =  f(t,p^).  (5.22) 

It  follows  that  g  depends  on  q,  ;,qp  only  through  the  value  of  (q-p): 

^  ^P2-Pi^   ^ 

g(t,q^,q2.P^.P2)  =  g(t,  q- P.P^^Pg)  =  g   t,q ^'Pi^P2  (5.23) 

V    IP2-P1I      y 

and  from  (5. 21)  we  may  easily  derive  the  identity 

~  I   -   I 

^   j^  g(t,q-P,p^,Pg)d(q.p)  =    ^^   f(t,p^)f(t,P2)  +  O(j^).       (5.24) 
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After  replacing  ff  by  g  in  (5.I6),  we  see  from  (5.I8)  that  we  may- 
replace  the  solution  operators  by  their  limiting  values  (for  t  -*  ") 
except  for  q^  in  the  spheres  A  and  B  (where  incomplete  collisions  occur). 
The  fact  that  the  cylinder  increases  in  length  with  t  (and  in  fact  be- 
comes infinite  for  t  -»■  "»)  is  rendered  harmless  by  the  fact  that  g 
vanishes  for  |q-p|  >  o"  |Pp-  P-,  I  +  ^-1  •   Thus  we  obtain 

L^  22  J 

The  term  0( — )  represents  the  error  due  to  incomplete  collisions. 

Although  we  have  derived  (5.25)  and  (5.2U)  by  an  argument  which  is 
valid  only  for  the  case  of  classical  mechanics,  we  shall  nevertheless 
proceed  to  show  that  the  Uehling-Uhlenbeck  form  of  the  quantum-mechanical 
Boltzmann  equation  follows  from  (5-2^)  and  (5-25)  provided  the  solution 
operators  and  the  operator  <?_  in  the  latter  are  given  their  quantum- 
mechanical  meaning. 

Let  G  be  the  operator  whose  Wigner-representer  is  g,  i.e. 

W(G)  =  g(t,q^,q2,P^.P2)-  (5-26) 

Then  from  (3-l6) 

W(Q^G)  =  l^g.  (5-27) 


Let 


a(t,q^,q2,Pi,P2)  =  lim  s[^'°^ s[%[^'''\e-^^B  =   W(a)  .     (5-28) 
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Then  from  our  definition  of  the   solution  operators    (3.12), 

iliT  ilnilw  ilw  ilff  iln 

2ft  ^2, o      "   ft     "2      2ft  "2, o   ^  ^      "    2h  ^2, o      ft       2  ^      2ft  "2,  o        ^  ^ 
A  =   lim  e  '      e  e  '      Q^G  e  '      e  e  '      -   Q^G- 

T->-oo 

(5-29) 


Here  H„  and  H„   are  the  Harailtonians  corresponding  to  s\        and  S.    ' 

d  d,0  t  T^ 


1,2    2, 


ki-^2l 


1  /  2 


2,0 


^2  =  a^J^Pi  -*-  P2) 


(5.50) 
(5.51) 


In  order  to  apply  the  results  of  scattering  theory  to  (5-29)  we  intro- 
duce center-of-mass  and  relative  coordinates: 


j(q^  +  qg)  ,  q  =  q2  -  ^1  ^  <ll 


^  >    ^2  "  '^c  "^  2  ^ 


Pc  =  Pi  ^  ^2  '  P  =  2^P2  -  Pi)  '  P 

M  =  2m  ,  \x  = 


,  =  i  P,  -  P  ,  P2  »  ^  P,  +  P  >     (5-32) 


2,H^  =  ip^E=E  +  «|,|),E  =  ^p2 


It  follows  that 


2,0 


H2  =  H  +  H^    ;    H^  ^  =  E  +  H^ 
From  (5.29), 


(5.53) 


lim 


(5.5^) 


where 


Here 
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TT  ilw  11  -a       2:1  -a 

2n  ^2,0      "    2fi  ^2     "   2fi  "2     2ii  ^2,0 


N(t)   =   e"^     ^'^  e     ^''     "e     ^''     ^e^'     "'"  =  M(t)M  (-t)  .      (5.55) 


2ft  "2,0      "    2fi  ^2 


M(t)    =   e  '      e  .  (5.56) 

Te^H  1    =  0  and   Ph^H  1    =  0,    it  follovs   that 


i^(E+H    )      -  ^(H  +H)  ±1  E      -  ^  H 

M(t)   =  e  e  =  e  e  .  (5-37) 


From  the  theory  of  time -dependent   quantiam-mechanical  scattering 


4,6,21 


22 
under  appropriate  conditions   on  the  'perturbing  operator'  (j)(|q|),  we  have 

the  following  result: 

lim  M(t)  =  SU  ;   lim  M(-t)  =  U.  (5-58) 

T  ->  00  T  -^  00 

Here  S  is  the  quantum-mechanical  scattering  operator  and  U  is  a  \xnitary 
operator  which  transforms  the  'unperturbed  Hamiltonian'  E  into  the  'per- 
turbed Hamiltonian' ,  H: 

H  =  U*EU   ;   U  U*  =  U*U  =  1.  (5-39) 

The  scattering  operator  is  also  unitary,  i.e. 

S  S*  =  S*S  =  1.  (5-^0) 

From  (5-38) 

lim  M*(t)  =  U^'s*  ;   lim  M*(-t)  =  U*.  (5-^1) 
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Thus  from  (5.55) 

lira  N(t)  =  S  U  U*  =  S,  (5.^2) 

T  -»  00 

and  from  (5.3^1-) 

A  =  S  Q^G  S*  -  Q^G.  (5.^3) 

It  is  somewhat  more  convenient  to  work  with  the  operator  J  =  S  -  1 
instead  of  the  scattering  operator  S.   Thus 

S  =  1  +  J,  (5-^^) 

and  from  (5 A3), 

A  =  J  Q^G  +  Q^G  /  +  J  Q^G  J*.  (5-^5) 

From  its  definition  (5. 21)  we  see  that  g  is  a  symmetric  function,  i.e. 
g(t,x  ,x  )  =  g(t,Xp,x  ).   From  this  it  follows  that  G  =  G(l,2)  is  a  sym- 
metric operator,  and  from  (5.^2)  and  (5.35)  we  see  that  S  and  J  are 
symmetric  operators,  since  E     and  Hp   are.   It  follows  from  (2.27)  that 
G  and  J  commute  with  Q^.   Since  Q  =  Q  and  Q  =  Q  ,  (5.^5)  becomes 

A  =  J  QgG  +  G(J  Q^)""  +  J  Q2(J  Q.^)* 

=  T  G  +  G  T*  +  T  G  T*,  (5-^6) 

where 

T  =  J  Q2.  (5.^7) 

From  (5.^0)  and  (5-^^) 

J  +  J*  =  -J  j"*^  =  -J*J,  (5.^8) 
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hence  from  (5-^7) 

T  +  T  =  -T  T  =  -T  T.  {^.^9) 

The  operator  J  is  most  conveniently  discussed  by  means  of  the  repre- 
senation  in  terms  of  a  complete  set  of  commuting  observalbes  that  includes 
the  unperturbed  Hamiltonian^  E;  hence  we  shall  introduce  that  representa- 
tion.  To  do  this,  let 

r  =  IpI  =  i  IP2-P1I   ;   P  =  r-\   =  -^^  .  (5.50) 

IP2-P1I 


Then 

P 

is 

a 

unit 

vector. 

and 

E 

2 
r 

\v,-v,f 

«,, 

(5.51) 

8^ 

From  (5-32)  one  may  determine  that  the  Jacobian  of  the  transformation  from 
p^^Pg  to  p,p^  is  1,  i.e. 

dp^dpg  =  dpdp^  (5.52) 

and  from  (5 '50)  we  see  that 

dp  =  r^drdA(p).  (5-55) 

Here  dA(p)  is  the  area  element  on  the  unit  sphere  |pj  =1.   From  (5. 55)  we 
obtain 

dp  =  X^dEdA(p)  (5-5^) 


where 


X  =  (2^5e)^/^  =  (nr)^/2.  (5.55) 


In  Section  2  we  introduced  the  momentiim  representation  R.   Let 
<^^i^>P,V^,Vy'-,vJ   =   X({)(p^,P2,p^,  ...,p^)  =   XR{f).  (5.56) 

Then  from  (2.8),  (5-52),  and  (5-51) 

2 
XR(E$)  =  ^|l  (t)(Pi^---^Pn)  =  ^'^^i^>P^V^,Vy-,V^),  (5.57) 

and  from  (2.5I),  (5-52),  (5-5^),  and  (5.56) 

II  I  f   =  r  |(t)|VdEdA(p)dp^dp^...dp^  =j"|(l)^|^dEdA(p)dp^dp^...dp^.   (5.53) 

(The  integrations  indicated  are,  with  respect  to  E  over  the  real  semi-axis 
0  S  E,  with  respect  to  p  over  the  unit  sphere,  and  with  respect  to  the  vec- 
tors p  ,p  , ...,p  over  all  space.)   It  follows  from  (5-56  -  5.58)  that  the 
correspondence  (5.56)  between  elements  ^  of  ^  and  functions 
^   (E,P,p  ,P,^..'^P  )  is  the  representation  in  terras  of  the  unperturbed 
Harailtonian  operator  E  and  the  operators  p,p  ,p  ,...,p  .   For  brevity  we 
denote  this  representation  by  R  .   Thus 

\(I)    -   '^J^,P,V^>Vy',V^)    =  XR(  $  ).  (5.59) 


Now,  in  terms  of  the  representation  R  ,  the  operator  J  is  given 


by 


^(J  $)  =   r  j(E,P,P')(t'o(E.P',P^.P3.---,P^)dA(p'),  (5.60) 


and  is  related  to  the  quantum-mechanical  differential  scattering  cross- 
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section,  ^(E,p,p'),  by 


21   1  |2 

Ti(E,P,p')  =  (2jt«)   r-^j(E,P,p')   .  (5.61) 


The  operator  J  can  be  put  in  the  form  of  an  integral  operator, 


Ro(j$)  =j^j(E,P,P')&(E-E')6(p^-p;^)6(p^-p^)...5(p^-p;)(|)^(E',p',p',p^,  ...,p;^) 


X  dE'dJV(p')dp^dp^...dp^,  (5-62) 


with  a  'singular  kernel'.   In  earlier  treatments,  representations  of  J  of 
the  form  (5.62)  were  used,  and  this  led  to  inconsistencies  involving  the 
square  of  the  6-function.   To  avoid  these  difficulties  we  shall  make  no 
use  of  (5-62)  and  shall  work  directly  from  (5-60). 

In  the  two  cases  Q  =  ^  ,    and  Q^  =  ^  ,  we  have,  from  (2.^4-2), 


ig  =  |(1  +  P^  g)  (5.65) 


where  P    is  the  permutation  that  interchanges  1  and  2.   From  (2.17) 

and  from  (5.56)  and  (5-50) 

Ro(^2^)  =  l[^o^^^P'Pc'P3'---'Pn)  ±(l)„(E,-P,p^,P5,...,p^r|.         (5.65) 
Thus,  from  (5-^7)  and  (5-60) 


k2 


=  I  y  j(E,P,p')^q(E,pSp^,P^,...,P^)  +  (t)^(E,-PSp^,p^,...,p^)JdA(p') 
=  y'b(E,P,p')(l)^(E,pSp^,p^,...,p^)dA(p'),  (5.66) 

where 

b(E,P,p')  =  I  h(E,P,P')  +  j(E,p,-p')J.  (5.67) 

Corresponding  to  these  two  cases,  we  define  the  symmetrized,  and  the  anti- 
symmetrized  cross-sections 

,2 

(5.68) 


a(E,P,P')  =  (2rt«^)|r'^b(E,P,p')|^  =  (2rtft)^U:[j(E,  P,  p' )+j(E,  P, -P' )]  I  . 


69) 


Of  course,  in  the  case  of  Maxwell-Boltzmann  statistics,  Qp  =  1,  T  =  J, 
b  =  j,  and  a  =  t)  . 
From  (5-66) 

R^(T*  $)  =  jrb''(E,P',p)(t)^(E,p',p^,Pj,...,p^)dA(p>),         (5. 
and  from  (5.^9)  we  obtain 


b(E,P,P')  +b  (E,P',P)  =  -  /  b(E,P,p")b  (E, P' , p" )dA(p" ) .      (5-70) 


For  p  =  p'  this  becomes 


b(E,P,p)  +b*(E,P,p)  =  -  /  |b(E,P,p')|^dA(p').  (5-71) 
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Our  object  now  is  to  express  the  operator  A,  as  given  by  (5-^6) 
first  in  R  -representation  and  then  in  R-representation.   From  the 
latter  it  will  be  easy  to  transform  to  Wigner  representation  and  then 
to  evaluate  the  right  side  of  (5.25)  using  (5.28).   In  order  to  do  this, 
let 

R(G)  =  7(p^,P2.p[,Pp  (5.72) 

and 

R^(G)  =  r^(E,p,p^,ESp;pp,  (5.73) 


\(^   $)  =/7'^(E,P,p^,ESPSp;,)(j)^(E',p',,p',Pj,...,p^)dE'dA(p')dp'.  (5.7^) 

Then,  from  (5-59),  (5-56),  (5-52),  (5-5^),  and  (5-72) 

XR(G  $)  =  Rq(G  f )  =   /  7'^(E,P,p^,E',pSp^)X'({)(p|,P^,P^,  ...,p^)(X')"^dp|dp^ 

=  X  /  7(v2_'V2''P{''p2^H'P[,V2'Vy-,V^)^V['^V'2'  (5-75) 

From  this  it  follows  that 


(XX')"S^(E,P,p^,E',pSp^)  =  7(p^,P2,p{.Pp 


=  r(|  p^-rp,  I  p^+rp,  I  p;-r'p',  |  p'  +  r'P')-   (5-76) 

Although  7   and  7  were  defined  as  the  R-  and  R  -  representers  of  G(1,2), 
it  is  clear  that  (5-76)  is  an  identity  which  holds  for  the  R-  and  R^- 
representers  of  an  arbitrary  operator,  B(1,2). 


Now  from  (5-66)  and  (5-73) 
R^(T(^)  =  /  b(E,p,p")7Q(E,p",p^,E',pSp^)(t)^(ESpSp^,p^,  ...,p^) 

X  dE'dA(p')dp^dA(p").  (5-77) 

Similarly  from  (5-69)  and  (5-75) 


R^(GT  ^)  =  /  b  (E',pSp")7q(E,p,p^,ESp",p^)(|)^(E',pSp^,P^,...,P^) 


X  dE'dA(p')dp'dA(p"),  (5-78) 


and 


R^(TGT*|)  =  /  b(E,p,p"')b*(ESpSp")7Q(E,p"',p^,ESp",p^) 


X  (j)^(E',p',p^,p^,  ...,p^)dE'dA(p')dp^dA(p")dA(p'").    (5-79) 


Let  R  (A)  =  a  (E, p,p  ,E',p',p').   Then  from  (5-^6); 


a^(E,p,p^,E',p',p^)  =  R^(A)  =  /  b(E,p,p")7^(E,p",p^,E',pSp^)dA(p") 


•  /  b*(E',p',p")7^(E,p,p^,E',p",p^)dA(p") 

•  Fid{E,p,p"'   )b*(ESpSp") 


7^(E,p'",p^,E',p",p^)dA(p")dA(p'").    (5.80) 


U5 

If  we  set  R(A)  =  a(p, ,p„,p' ,p' ),  then  (5-76)  is  valid  for  a  and  a  as 
J.      d.      X      d.  o 

well  as  for  7   and  7    .      Thus  from  (5.76)  and  (5.80) 
a.(Pj_,P2.P{.P2)  =  R(A)  =  /  b(E,  p,  p")7(2P^-rp",ip^+rp",p^,ppdA(p") 

+  /  b  (E',pSp")7(Pj^,P2.2P^-r'p",?p^+r'p")dA(p") 

+  /  b(E,p,p"')b*(E',p',p") 


X  7(?p^-rp'"  ,^^+rp"Sip^-r  •  p",^^+r '  p"  ) 


XdA(p")dA(p'").  (5.81) 

By  definition  7  =   R(G)  and  g  =  W(g) .   Hence  from  (5-3) 

,       .   .^    r  sfrl*^Pi"Pl^-*-^2-(P2-P2^]  r        ,.     ,.  ,.     ...    . 
7(p^,P2.P-|_.P2)  =  /  e  gLt,s^,S2,2(p^+p|),2(P2+P2)J'is^c 

(5.82) 

We  introduce  new  variables  of  integration,  s  =  Sp-  s  ,  s  =  -^(s  +  Sp) . 

Then  ds^ds^  =  dsds  ,  and 
12       c' 

,        ,   ,^    r  ^K'^P^P2-VP2)-^*^-(Pl-Pi^P2-P20 

X  g[i(p^+p|),i(p2+P2)]<isds^.  (5.85) 
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Here  we  do  not  explicitly  Indicate  the  dependence  of  g  on  t,  s,  and  s  . 
But  from  (5.23)  g  is  independent  of  s  ;  therefore 

X  g\j{Vj_+v{),kv2+v'2)]ds.  (5.84) 

If  we  insert  (5-84)  in  (5.81)  we  obtain 
a(p^,P2,p{,pp  =  (2nfi)^6(p|+p^-p^-P2)  i  J  lD(E,p,p")e*^ 

X  g[j(2p|+r'p'-rp"),i(2p^-r'p'+rp")JdsdA(p") 


-  s.(r'p'-rp  ) 


+  /  b*(ESp',p")e* 


^  s.(r'p"-rp) 


:  gJi(2p^+rp-r'p"),i(2p2-rfH-r'p")JdsdA(p") 

i    /  1  't    t.i\ 
p  ^  T-  s.(r'p  -rp"') 

•  /  b(E,p,p'")b  (ESpSp")e 


X  g[i(p^-rp'"-r'p"),i(p^+rp"'+r'p")]dsdA(p")dA(p"')i  .(5.85) 

By  definition  a  =  R(a)  and  a  =  W(a) .   Therefore  from  (5-2) 
-2i, 


-6  r  ^'^w^2-^-2) 


i(t,q^,q2,p^,P2)  =  (rth) 

X  ci{^^-z^,^^-z^,Vj+'Z'-^,V2+^2^^^1^^2'  (5.86) 


k-J 


and  from  (5-25)  and  (5-28) 


D,f(t,p,)  =  ^  W-^^f 


1  .  .^-6  r  V(^i-vi2-"2) 


X  <x{p^-z^,-p^-z^,p^+z^,j,^+z^)dz^d.z^dq^d.-p^  +   0(:^).         (5-87) 


The  integrations  with  respect  to  qp  and  z  may  be  carried  out  in  (5.87)' 
The  result  is 


-2i 


.f(t,p^)  =  ^  Ufi)-^Je  ^        1  a(p^-z,p2,p^+z,p2)dzdp2  +  0(:^) .    (5-88) 


We  now  insert  (5.85)  in  (5-88).   In  so  doing  we  replace  p  by  p,-z, 

P-[  by  p-|_+z,  and  p^  by  p^.   Then  5(pj^+p^-p^-P2)  becomes  S(2z)  =  2"^6(z). 

The  result  is 


s.(p-p') 


X  dsdA(p')dp2  +  /  b  (E,p,  p')e 


t,s.i£tP:i,p+|(p-p>),P2-|(p-P'^ 
|f>»-p'|  J 


(P-P-) 


X  g 


;,s.i£i:£:l,P^+|(p-p'),P2-|(p-p' 
|f>t-p'| 


dsdA(p'  )dp^ 


■/ 


s-(p'-p") 


+  /  b(E,p,p")b  (E,p,p')e 


t,s-^P'-^P"^^-|(p"4-p'),^  +  j(p"+p'^dsdA(p-)dA(p")dpi  +  0{^)^ 
1p'+p"I  J  J 

(5.89) 
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In  (5-89)  we  have  indicated  explicitly  the  dependence  of  g  on  t  and  s. 
This  dependence  is  obtained  from  (5.23)- 

From  (5.32),  Pg  =  p  +  2p.   Therefore,  for  fixed  p  , 

dpg  =  2^dp  =  2^r^drdA(p).  (5-90) 

If  we  use  (5-90)  to  replace  dp_  in  (5.89),  we  may  then  use  the  result 
(A.1),(A.36)  of  the  Appendix  to  simplify  (5.89).   This  yields 


D^f(t,p-L)  =-^^^y  r"2b(E,P,p)g|t,s.p,p^,P2]d(s.p 
+  /  r"^*(E,P,p)grt,s-p,p^,p2d(s-p)dp2 
+  /  r"T)(E,P,p')b*(E,P,p')g[t,s-p',5^,P2] 


)dp2 


X  d(s-p')dA(p')dp^>+0(:i^).  (5.91) 


^    p  ^    p 

P^  =  ^  -  rp'  =  p^+  r(p-p'),   pg  =  -^  +  rp'  =  p^-  r(p-p').        (5-92) 


Let  p  =  pCpp-p-,)  =  rp".   From  (5.92)  we  easily  determine  that 

/S.V  /n2/^222 

P-]_  +  Pp  =  ?■[_  +  Pg  a-^<i  ?!     +  P2     =  ?!     +  P2    •      '^^^   ^1   ^'^  P2   ^^^   ^^'^ 

final  momenta  after  a  collision  of  two  particles  whose  initial  momenta 


direction  is  p'.   The  collision  preserves  total  momentum  and  the  kinetic 
energy  of  both  the  total  and  relative  motions. 


^9 


We  now  use  (5-71)  to  simplify  (5.91).   Thus  we  obtain 
D^f(t,p^)  =i2^J'r-2|b(E,p,p')|2M'g(t,s.p',p^,P2)d(s.p>) 

-yg(t,s-p,Pj^,P2)d(s-p)l  dA(p')dp2  +  0(i).  (5.95) 

Finally,  from  {'?.2k)    and  (5-68), 
D^f(t,p^)  =i  y'o(E,P,p')[f(t,9^)f(t,P2)-f(t,p^)f(t,P2r] 


'2   ^1 


dA(p')dp2  +  0{^)^ 


(5.9^^) 


For  sufficiently  large  1,  we  may  neglect  the  remainder  term  in  (5-9^)  and 
the  result  is  the  Uehling-Uhlenbeck  form  of  the  quantiMi-mechanlcal  Boltz- 
mann  equation,  for  the  spatially  homogeneous  case.  The  Uehling-Uhlenbeck 
equation  for  the  inhomogeneous  case  is  given  by 

D/(t,q^,p^)  =  i   ra(E,p,p')[f(t,q^,5i)f(t,q^,P2)-f(t,q^,p^)f(t,q^,P2)] 


'2   ^1 


dA(p')dp. 


(5.95) 


The  derivation  of  (5-95)  would  require  further  work. 


Conclusion: 


It  is  quite  clear  that  the  weak  point  in  the  present  derivation  is  the 
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argTiment  leading  to  (5-25)  and  {'^.2h)   which  proceeds  by  analogy  with  the 
classical  case.   In  fact,  in  our  opinion,  this  point  is  so  weak  that  it 
leaves  the  validity  of  the  result,  (5 -9^);!  in  considerable  doubt.   In 
recent  years,  some  work  has  been  done   to  clarify  the  asymptotic  approach, 
(as  t  -»  00  )  of  the  operator  exp(-^  E)exp(— ^  H)  to  the  operator  SU,  as  in 
(5-38).   It  may  be  that  this  type  of  study  will  resolve  the  difficulty 
that  we  have  pointed  out. 
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Appendix 

The  purpose  of  this  section  is  to  apply  the  formula  for  the  inver- 
sion of  Fourier  transforms  in  order  to  simplify  integrals  of  the  form 

ir 


'/=" 


;.(p-w)  r- 


P+w| 


dA(p)dA(co)ds. 


(A.l) 


Here  p  and  co  are  unit  vectors  and  dA  is  the  area  element  on  the  unit 
sphere.   We  shall  make  use  of  the  Fourier  inverse  relation  in  the  form 


/ 


+—  y(x-x') 
;   ^        f(x')dx'dy  =  ^  f(x), 


(A.2) 


The  area  element  dA  can  be  replaced  by  an  area  element  in  the  plane 


by  introducing  a  cartesian  coordinate  system  such  that  p  has  components 

,222 
WPj^'  P^4-P2+  P3 


2    2    2 
p  =  (p  ,p  ,p  );  p  +  p  +  P  =  1.   Then 


dA(p) 


^1 


(A.3) 


To  transform  the  integral  I,  we  introduce  a  vector 


T  =  p  -  U) 


(aA) 


and  a  unit  vector 


|p  +  w| 


(A.5) 
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Then 

a-T  =  0.  (A. 6) 


Let 


3^  =   (a^+  al)-^^^{o^,-a^,0),  (A.7) 


/  2    2>-l/2,  2  2-  ,^  OS 


jg  =  a.  (A. 9) 

j,  ,  j  ,  and  j^  form  an  orthonormal  basis.   Since  o-x  =  0, 


T    =    Ti^jj^    +    T^j^  (A. 10) 


where 


Ji^  =  B"^/2[a2(Pi-u)^)    -    a^(P2-a)2)],  (A-ll) 

^^^[v3^Vl)-*-''2''3^V2)-^^V3[|'  (A. 12) 


Let 


X5   =    T.j^   =   B 


B  =   al+  o\.  (A. 15) 


A  =    (p+w)^  =   2(l+p-co)    =   2(l+p^a)^+P2   g+P  u)^) .  (A.l4) 


Then 


j^  =  A"l/2(p^-HJ^),  ag  =  A"^/^(P24^2)  ^^-^5) 


-iF  2  2      2  2I 

B  =  A        p^+2p^a)-|^+u)^+P2+2p2CJ2+W2    " 


(A. 16) 
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We  proceed  to  calculate  the  Jacobian  of  the  transformation  from 
p  and  u)  to  T  and  o,  evaluating  the  results  only  at  t  =  0,  i.e.  at 
p  =  cj  =  g.   By  using  relations  of  the  form 


ap^  P5 

we  successively  obtain,  at  t  =  0, 


(A. 17) 


A  =  h, 

hA 


hp 


hp. 


SP, 


^^1. 


SP. 


=   B 


[1 

3 

^1 

-1/2^ 

■1/2, 


SP. 


Sco-, 


=  0, 


Bp 


Bt, 


ap. 


St, 


Co,, 


^T. 


SP. 


Co. 


2     ^'"2 

B-1/2. 
-B    a^, 

.-1/2. 


Sw^ 


•Ca, 


St, 


Su). 


(A. 18) 


(A. 19) 


(A. 20) 


(A. 21) 


(A. 22) 


(A. 25) 


(A.2i.) 
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•1/2. 


a/2   -1 


(A. 25) 


da-j^dagd-Ti^dT^ 


5a,   5a,   5o,   5a, 

Sag  Sog  Sag  Sag 
5p^  Spg  So)^  Swg 


^4  ^^4 


^\     ^\ 


5p  Sp  Soj  5w 
5t  St  St  St 
Sp^     Spg     S(i)-j^     Su)g 


(A. 26) 


If  ve  insert   (A. 20),    (A. 21),    (A. 25),    and  (A.2J+)    in   (A. 26)   we   obtain  easily, 
at  T  =   0, 


(A. 27) 


Let  J  be  defined  by 
dA(p)dA(co)   =  JdTj^dT   dA(a) 
Then  from  (A. 3), 


(A. 28) 


dp^dpgdto^dcog                           do-j^dog 
=    J    dTj^dT^ 


IP5I     hjl 


(A. 29) 
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and,  at  T  =  0, 


l°3l       -1  .     ^ 

J  =  2 J  -^  =  1.  (A. 50) 


IP5I  I-3I  ' 


Now  define  F  to  be  zero  in  (A.l)  for  |p|  ?^  1  or  |u)|  5/  1.   The  inte- 
gral may  be  written  as  a  sum  of  four  terms  depending  on  the  domain  of 
integration, 

I,  fp^  >  0,co,  >  0),I   (P,  <  0,u)^  <  0),I   (P,  >  0,0)^  <  0),I   fp,  <  0,cj,  >  0) 

(A. 51) 

and  by  (A. 28) 

ir 
I  =  /  e^    F[s-a,P,u)JjdTi^dT  dA(a)ds.  (A. 52) 

In  the  transformed  terms,  I   corresponds  to  a,  >  0,  and  I   to  o-,  <  0, 
++  5  —     0 

and  the  T-domains  of  both  of  these  integrals  include  the  point  t  =  0. 
The  T-domains  for  I  and  I  exclude  the  point  t  =  0-  (A. 52)  may  be 
written 


p  -T-fs.  T,  +s  T  1 
I  =  /  e^L^    5  ^JF[sg,P,a)]jdSgdSj^ds  dT^dT  dA(a),  (A. 55) 


and  from  (A.2) 


I  =  ( 


%^)7^[v^^^]^ 


dA(a).  (A. 5^) 

T=0 
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But  atT  =  0,  p  =  u)  =  a  and  J  =  1;  hence 

I  =  i^)^     /  Frsg,a,oldSgdA(o).  (A.  55) 

Since  the  point  t  =  0  is  outside  the  domain  of  integration,  ^4.  =  -'-  ."  *^* 
The  contributions  I   and  I   may  be  added  to  produce  an  integral  over 
the  entire  unit  sphere  |a|  =  1.   By  changing  o  to  p,  we  obtain  the  result 

o 

1  =   i^)         r  F[(s.p),p,p]d(s.p)dA(p).  (A. 56) 
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